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Abstract 

Hysteresis has been observed in a study of the transition between laminar 
flow and vortex shedding in a quasi-two dimensional system. The system is 
a vertical, rapidly flowing soap film which is penetrated by a rod oriented 
perpendicular to the film plane. Our experiments show that the transition 
from laminar flow to a periodic Karman vortex street can be hysteretic, i.e. 
vortices can survive at velocities lower than the velocity needed to generate 
them. 
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Pattern forming fluid flow is often difficult to understand even at quite modest velocities. 
A widely studied example is the chain of vortices generated downstream from a rod placed 
in a uniform flow, with the axis of the rod oriented perpendicular to the flow direction. In 
a pioneering study of this problem, Karman analyzed the stability of this vortex street [|TJ . 
His starting point was the assumption that the flow is two dimensional, i.e., no velocity 
variations are generated along the rod axis. He was concerned only with the stability of 
the vortex array and made no attempt to explore the transition itself. Although researchers 
have put considerable effort to better understand this phenomenon, which is considered to 
be a precursor to turbulence, a full understanding is still lacking 0. 

It is now widely accepted that the transition (primary instability) from laminar flow 
to vortex shedding is supercritical in nature; on changing the mean flow speed V around 
the transition point, vortex shedding appears or disappears at a unique velocity V c . Below 
(above) V c the only stable state of the system is that of the laminar (vortex shedding) flow. 

In an important wind tunnel experiment on vortex shedding, Provansal et al. || showed 
that the primary instability occurs at a critical Reynolds number Re c ~46. The vortex 
shedding disappeared at the same velocity at which the onset occurred, i.e. hysteresis was 
absent. These measurements are well described by the Landau-Stuart theory for supercritical 
transitions. Sreenivasan et al. @ observed that an acoustic disturbance or a vibration of 
the rod itself could activate the vortex street when Re was less than Re c but greater than 
some minimum value Re m in — Re c /2. Their experiment showed that the general behavior 
of the system was position independent. The Landau-Stuart theory was adequate to explain 
their experimental findings too. Couder || and Gharib || were among the first to study 
vortex shedding behind a cylinder in soap films. In particular, Gharib and his associates 
carried out systematic investigations of the vortex shedding frequency / as a function of V 
in a horizontally flowing film 0. 

In this paper we present measurements demonstrating that the primary instability can 
be hysteretic. This finding is also in contrast with numerical results J7J. Although hysteretic 
transitions are well known in fluid dynamics ||, we know of no published laboratory ex- 
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periments or computer simulations, which suggest that the onset of vortex shedding from a 
rigid rod is a hysteretic phenomenon. 

Our experimental setup consists of a rapidly flowing soap film formed between two verti- 
cally positioned 0.25 mm diameter nylon lines separated by 6 cm (see Fig. 1). The lines are 
parallel along a 45 cm vertical segment. They are joined at both ends of the segment. The 
film is fed continuously from the top through a high precision metering valve. The valve 
was opened and closed at different rates by a motor connected to it, thus the acceleration 
dV/dt was under experimental control. The film is a mixture of soap (Dawn) and distilled 
water with a typical soap concentration of 1 wt.%. 

Using an infrared absorption technique, the soap film thickness profiles parallel and 
perpendicular to V were measured. In the absence of the rod, the horizontal film thickness 
was essentially flat with less than 5% variations, whereas there was systematic thinning 
along the vertical direction. Due to a small acceleration of the film, the latter amounts to 
a 1% change in the thickness over a distance between the two LDV probes. The average 
thickness, which was found to be several microns, does not change over the range of V 
spanned by our experiments. 

The velocity measurements were made using a dual head laser Doppler velocimeter 
(LDV). At the center of the parallel segment of the soap film, a glass rod of diameter d=l 
mm penetrates the film in the z direction. The velocities V and the fluctuating horizontal 
component V x (t) were measured at points 5 cm above and 5 cm below the rod respectively. 

At a certain critical flow rate there is a transition from the laminar flow (LF) state to 
the vortex shedding (VS) state. Slightly below this critical velocity, the flow is characterized 
by two counter- rotating vortices underneath the rod (see Fig. 2a). In the VS state these 
recirculating vortices peel off the rod and flow downstream. The continuously generated 
counter-rotating vortices form a vortex street (see Fig. 2b), as described by Karman [![]. 

All measurements were made near the transition. The experiments start with a value 
of V corresponding to the LF regime. The valve is then opened to slowly increase V. The 
critical velocity where VS commences is called V^ p . After some waiting time the valve 
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is slowly closed. As a result V(t) decreases, and at V{t)=V down the system undergoes a 
transition from VS into LF. This cycle was repeated several times in each run. To find the 
transition from LF to VS, V x (t) was recorded together with V. In theory V x (t) is constant 
in the laminar state and equal to zero directly below the center of the rod. In the VS regime 
V x (t) is expected to display periodic behavior. 

In order to determine the transition velocity accurately we have calculated the velocity 
probability distribution function P(V X ). For this purpose the total recording time was 
segmented into intervals At ~200ms. Then P(V X ) was calculated within each interval. 
The intervals At were much smaller than the characteristic changing time of V(t) but large 
enough to include a statistically significant number of data points and many periods of 
oscillation in the VS state. A typical result is shown in Fig. 3. Here the magnitude of 
P(V x ,t) is mapped into gray scale values and shown as a function of V x and t. It can be 
seen that sharp changes in P(V X , t) precisely indicate the transition between the LF and VS 
regimes. In the LF regime, P is sharply peaked at V x =0, but in the VS regime P has a 
broad distribution, as for a (noisy) sine wave. These measurements provide a powerful way 
to accurately determine the times t* d and t* of LF— >VS and VS— >LF transitions ||10|| . 

The hysteresis is apparent in Fig. 3, where V(t) and P(V x ,t) are displayed together. 
By determining the value of V(t) at the transition times t* and t* d one find that V^ p = 
V(t* u ) is not equal to V c down = V(t* d ), that is, there is no unique critical velocity for the 
transition. In this analysis we have incorporated the fact that V is noisy at the level of 
±4%. Therefore to determine V down (and V^ p ) we always used the lowest (and highest) 
values of (the fluctuating) V at the transition. 

The critical velocities were measured in 20 experimental runs. The ensemble average of 
the different experimental runs provide < V c down >= 0.42 ± 0.02m/s. The relative difference 
(< V™ p >-< V down >)/ < V down > is approximately 14%. The Reynolds number -which is 
defined as Re = Vd/v- is roughly 50 in the hysteretic gap. A large uncertainty of this value 
comes from the poorly determined two-dimensional soap film viscosity u, which depends on 
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the film thickness HIT . 



Next, we demonstrate that the observed hysteresis is not a result of some delay in the 
response of the system to the bifurcation, as it would be the case for dynamical hysteresis 
I2| , |T3|| . We set the mean flow rate constant in the interval [V^, n ,V^* p ]. In this case the 
laminar state can persist for an indefinite length of time. However applying sufficiently 
large acoustic or pulse-like mechanical perturbations, the system could be driven into the 
VS state. The system remains in this state even if the perturbation is turned off. Applying 
similar perturbations VS— >LF transition was never observed at constant V. As a result we 
can conclude that we observed static hysteresis, i.e. its existence is independent of the rate 
of change of V. 

At the transition, a certain time interval is required for the new steady state to establish 
itself as it is seen in the close vicinity of t* in Fig. 3. The overall time of the transient remains 
finite even if V is changed very slowly. This would not be true for continuous transitions 
||16|| . Careful examination of P(V x ,t) shows that the most probable V x grows exponentially 
at the transition points. Using exponential curve fitting one finds that r 9rowth ~ 1.32± 
0.05 s. Thus it was necessary to vary V slowly enough, that this finite transit time did not 
interfere with our identification of V^ p and V^ own . 

Because hysteresis is not seen at the onset of vortex shedding in 3D systems, we have 
looked for possible side effects that might be responsible for our unexpected finding. It is 
known that oscillations of the rod can induce hysteresis ||. Our rod is rigidly supported, 
so its vibration can only be induced by the fluid itself |J. Since the film is only a few 
microns thick, it carries insufficient momentum flux to deflect our 1 mm diameter glass rod 
more than a nanometer. The possible motion of the vertical fishing lines which support the 
film was also considered. Taking into account the geometry of our experimental setup, we 
estimate that their deflection, produced by the passage of an eddy, is no more than a micron. 
Auxiliary experiments were also performed to exclude wetting properties and the effect of 
the air boundary near the film as possible sources of hysteresis PO . 



Our experiments are the first to establish that the transition to vortex shedding can be 



subcritical. In the absence of any available theory, we consider a generic model based on 



an amplitude equation fn| , l5|] , namely a fifth order Landau equation [[HJ for the complex 
velocity v—v x + iv y : 

dv{t)/dt = av + 2bv 2 v — cv 4 v, (1) 

where all variables indicated with tilde are complex and Vi = Vi~ < Vi >t (i=x or y). 
The complex velocity can be written as v=ve lr ^, where 4> is the phase factor and v is the 
amplitude. In this case the evolution of v is governed by the real part of Eq. ([!]): 

dv/dt = v (a r + 2b r v 2 — c r v 4 ), (2) 

where the subscript r is used to designate the real part of complex numbers and a r is 
considered to be proportional to (V — V^ p ). In the experiments a r is adjusted by varying V 



In the following we consider the stationary solutions of Eq. (|2|) . The laminar solution is 
vlf = 0. Note that at the position where V x is measured, (V x ) t =0. Therefore v x (t) = V x {t) 
for that particular point. The time-dependent solution is described by: 
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1± Jl + a r c r /b 2 r ). (3) 



It suffices to retain only the positive values of Eq.(^). Since the minus sign in this equation 
corresponds to an unstable state, the observable vortex shedding is represented by Vy S . 

After some elementary algebra one can draw the stability (or bifurcation) diagram of this 
model (see Fig. 4, where only the positive solutions are shown). The solution vlf is stable 
only if a r is smaller than a critical value (a r )" p =0 (see regimes I and II). The solutions Vy S 
exist only for a r >(a r ) d c own , where (a r )f wn =-b 2 r /c r . 

If a r < (a r )c° wn then the system is in the LF state. On increasing a r , the system 
enters into the hysteretic regime (II), where the system can stay in the LF state as long 
as the amplitude of the fluctuations does not exceed the value of Vy S . If a r > (a r )^ p 
then the LF state becomes unstable; even infinitesimally small noise will force the system 



to switch from vlf to Vy S . In regime II, finite perturbations with amplitude larger than 
v vs force the system to switch to the VS state. Such a switch is clearly inconsistent with 
a supercritical bifurcation, for which the saturation amplitude continuously changes with 
the control parameter. The everywhere stable solution Vy S does not exist below (a r )f own . 
Therefore the system must return from the VS state to the LF state at a r = (a r )f> wn even 
when noise is absent. In this sense the VS^LF transition is more robust than its inverse. 

Our observations are consistent with all the above described features of the model. How- 
ever the application of this model for data fitting purposes is not obvious, since the rela- 
tionship between the parameters a r ,b r ,c r and the control parameter of the experiment is 
unknown. 

We use the above model and choose a r to be a r = c^e, where e = (V — V^ p )/V^ p . 
Other parameters are kept as fitting constants. In this case the relevant fitting parameters 
are ai/c r , b r /c r , and V c ap . Figure 5 shows the results of our fits together with the same 
experimental data set presented in Fig. 3. Shown here is the ^-dependence of the values of 
v x for which P(v x ) is maximum. This figure demonstrates both the sharpness of the VS— >LF 
transition and the breadth of the LF— >VS transition in our experiment. Open symbols show 
the measurements with decreasing V and closed symbols refer to measurements made in the 
opposite direction. 

The dotted line in Figure 5. is the best fit to the data with V c fixed at a value of 0.504 
m/s, which is the largest mean velocity where the laminar state was still observed. This fit 
seriously underestimates the width of the observed hysteretic gap. The dashed line is the 
best fit with a different constraint: it was required that the fit must terminate exactly at 
the lowest velocity for which the VS state exists. In this case the part of the curve which 
represents the unstable trajectory, terminates on the laminar branch in the middle of the 
hysteresis loop, predicting that LF state should not exist above V=QA5m/s. The solid line 
designates a fit where both of the above constraints are imposed. 

Although the generic model of the hysteresis is in fair agreement with our observations, 
none of the above fits are very satisfactory. It is interesting to note that a simple three- 
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parameter phenomenological equation Po(V — V c ) + 2P\V — P2V 2 = provides a surprisingly 
good fit to our data in the VS state. Using least square fitting to the solution of this equation 



provides the following result: v x = 0.04 ± y V^/39.06 — 0.0113, where the velocities are in 
units of m/s. 

In summary, we have observed hysteresis at the onset of the Karman vortex street in a 
quasi two-dimensional film. In view of prior studies of vortex shedding this result is sur- 
prising. Although the origin of the hysteresis is still not clear, we can exclude mechanical 
vibrations, wetting properties of the rod, and the air boundary layer as possible sources. 
Our preliminary experiments suggest that the effect may be connected with the fact that 
(a) the film is very thin in the recirculation region just before vortex shedding starts and 



(b) the kinematic viscosity of the film depends on its thickness [0,[Tj|,[nj . The fifth order 
amplitude equation is in qualitative agreement with our observations, but a simple phe- 
nomenological equation provides a better numerical fit to the experimental data in the VS 
regime. Theoretical explanation of our experimental findings as well as the full description 
of the vortex shedding continues to offer a great challenge for the future. 
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and by additional support from the Hungarian Science Foundation grant OTKA F17310. 



8 



REFERENCES 

[1] T. von Karman, Gott. Nachr. , 24 (1911). 

[2] There has been a vast amount of research on vortex shedding behind bluff bodies it is 
therefore not possible to review the literature in its entirety. For an excellent review see 
C. H. K. Williamson, Annu. Rev., Fluid Mech, 28, 477 (1996). 

[3] M. Provansal, C. Mathis, and L. Boyer, J. Fluid Mech. 182, 1 (1987). 

[4] K. R. Sreenivasan, P. J. Strykowski, and D. dinger, in ASME Forum on Unsteady 
Separation, edited by K. N. Ghia, volume 52, page 1, 1986, and private communication. 

[5] Y. Couder and C. Basdevant, J. Fluid Mech. 173, 225 (1986). 

[6] M. Gharib and P. Derango, Physica D 3, 406 (1989). 

[7] R. D. Henderson and D. Barkley, Phys. Fluids 8, 1683 (1996). 

[8] D. Brika and A. Laneville, J. Fluid Mech. 250, 481 (1993). 

[9] M. Rivera et al., Rev. Sci. Instr. 69, 3215 (1998). 

[10] In identifying the transition time it must be taken into account that V(t) and v x (t) is 
measured different places. Since V(t) near the transition is ~50 cm/s, the related time 
delay (about 0.2 s) is negligible to the characteristic changing time of V. 

[11] M. Rivera et al. (private communication). 

[12] G. H. Goldsztein, F. Broner, and S. H. Strogatz, J. Appl. Math. 57, 1163 (1997). 
[13] N. Berglund and H. Kunz, Phys. Rev. Lett. 78, 1691 (1997). 
[14] L. D. Landau, C.R. Acad. Sci. U.R.S.S. 44, 311 (1944). 
[15] J. T. Stuart, J. Fluid Mech. 4, 1 (1958). 

[16] S. H. Strogatz, Nonlinear Dynamics and Chaos, Addison Wesley, New York, 1994. 

9 



[17] The mean velocity V is used here rather than the Reynolds number because of an 
uncertainty concerning the value of the viscosity of the flowing soap film. 

[18] M. A. Rutgers et al, Phys. Fluids, 8, 2847 (1996). 

[19] P. Vorobieff and R.E. Ecke, Phys. Rev. E, 60, RC2953 (1999). 

[20] Z. Voros, W.I. Goldburg, V.K. Horvath, J.R. Cressman, and X.L. Wu, to be published. 



10 



Figure captions 

1. The schematic of the experimental setup. The diameter of the rod is 1mm. 

2. Interference pictures of the flowing soap film below the rod in (a) the laminar state 
and (b) the vortex shedding state. It was shown previously that the thickness of the 
film is correlated with the film velocity [10]. 

3. Here the function V(t) is shown together with the P(v x , t) map at |dV/<i£|=0.004 m/s 2 . 
The gray scale corresponds to the probability density of v x . The solid and dashed lines 
are intended to guide the eye for the transitions VS^LF and LF^VS. 

4. The bifurcation diagram of Eq. (fj). The turning point of the nontrivial solution is 
indicated by A. The solid and dashed lines designate the stable and unstable solutions 
respectively. 

5. The main figure shows our experimental data together with 3 different fits using Eq. 
(S). The circles and squares represent two different data sets. Closed and opened 
symbols are used for data taken at increasing and decreasing V respectively. The solid 
and dashed lines have the same meaning as in Figure 4. The inset shows our best 
parametric fit to the data (see text). 
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